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ABSTRACT 

Lya absorption spectra of QSOs at redshifts z ~ 6 show complete Gunn-Peterson 
absorption troughs (dark gaps) separated by tiny leaks. The dark gaps are from the 
intergalactic medium (IGM) where the density of neutral hydrogen are high enough 
to produce almost saturated absorptions, however, where the transmitted leaks come 
from is still unclear so far. We demonstrate that leaking can originate from the lowest 
density voids in the IGM as well as the ionized patches around ionizing sources using 
semi-analytical simulations. If leaks were produced in lowest density voids, the IGM 
might already be highly ionized, and the ionizing background should be almost uni- 
form; in contrast, if leaks come from ionized patches, the neutral fraction of IGM would 
be still high, and the ionizing background is significantly inhomogeneous. Therefore, 
the origin of leaking is crucial to determining the epoch of inhomogeneous-to-uniform 
transition of the the ionizing photon background. We show that the origin could be 
studied with the statistical features of leaks. Actually, Lya leaks can be well defined 
and described by the equivalent width W and the full width of half area Wh, both 
of which are less contaminated by instrumental resolution and noise. It is found that 
the distribution of W and Wh of Lya leaks are sensitive to the modeling of the ioniz- 
ing background. We consider four representative models: uniform ionizing background 
(model 0), the photoionization rate of neutral hydrogen Fhi and the density of IGM 
are either linearly correlated (model I), or anti-correlated (model II), and Fhi is cor- 
related with high density peaks containing ionizing sources (model III). Although all 
of these models can match to the mean and variance of the observed effective optical 
depth of the IGM at z ~ 6, their distribution of W and W-a are very different from 
each other. Consequently, the leak statistics provides an effective tool to probe the 
evolutionary history of reionization at z ~ 5 — 6.5. Similar statistics would also be 
applicable to the reionization of He II at z ~ 3 

Key words: cosmology: theory - intergalactic medium - large-scale structure of the 
universe 



1 INTRODUCTION 

In the last decade, the Lya forests of QSO's absorption spec- 
tra at redshifts z < 5 have played an important role in 
understanding the diffuse cosmic baryon gas and the UV 
ionizing photon background, and constraining cosmological 
models and parameters (e.g. Ranch et al. 1997; Croft et 
al. 2002; Bofton et al. 2005; Seljak et al. 2005; Jena et al. 
2005; Viel et al. 2006). Recently, more and more UV photon 
sources, including QSOs, GRB, Lyman-break galaxies, and 
Lya-emitters at redshifts z > 5 have been observed (see El- 
lis 2007 and reference therein). Due to the rapidly increase 
of Gunn-Peterson (GP) optical depth at z > 5, their absorp- 
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tion spectra show long dark gaps on scales of tens of Mpc 
separated by tiny transmitted leaks. It has been suggested 
that we are observing the end stage of reionization (Fan et 
al. 2006). 

It has been known that the dark gaps are from the IGM 
where the density of neutral hydrogen are high enough to 
produce almost complete absorptions, however, where the 
transmitted leaks come from is still unclear. In photoioniza- 
tion equilibrium, the density of neutral hydrogen nm oc ap^, 
here a is the recombination rate and p is the density of IGM; 
therefore, even when most of the IGM are neutral enough 
to produce complete Lya absorptions, it is still possible for 
the lowest density voids to provide prominent transmitted 
fluxes. On the other hand, the leaks can also come from ion- 
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ized patches around ionizing sources where the intensity of 
UV radiation are higher than average. 

The origin of leaking is crucial to understanding the 
history of reionization. According to commonly accepted 
scenario of the reionization, in the early stage, the ion- 
ized regions are isolated patches in the neutral hydrogen 
background (e.g., Ciardi et al. 2003; Sokasian et al. 2003; 
Mellema ot al. 2006; Gncdin 2006; Trac & Ccn 2006). and the 
subsequent growing and overlapping of the ionized patches 
load to the ending of reionization (e.g., Ciardi et al. 2003; 
Sokasian et al. 2003; Mellema et al. 2006; Gnedin 2004). 
Thus, if leaks mostly come from ionized patches, reioniza- 
tion should happen in the early stage. In contrast, if they 
were produced in lowest density voids, the the UV ionizing 
background might has already underwent an evolution from 
highly inhomogeneous to uniform distribution. 

A variety of statistics has been used to study the evolu- 
tion of reionization, such as the mean and dispersion of GP 
optical depth, the probability distribution function (PDF) of 
the flux, and the size of dark gaps (e.g.. Fan et al. 2002, 2006; 
Songaila & Cowic 2002; Paschos & Norman 2005; Kohlcr 
et al. 2007; Gallerani et al. 2006; Becker et al. 2006), but 
all of them seem to be ineffective to provide the informa- 
tion of leak's origin and the inhomogeneity of UV ionizing 
background. The GP optical depth is an average, and not 
sensitive to details of reionization. The statistical properties 
(mean and variance) of the GP optical depth at a ~ 6 can 
be well explained by either the fluctuation of ionizing back- 
ground (Fan et al. 2006) or models with uniform ionizing 
background (Lidz et al. 2006; Liu et al. 2006, hereafter Pa- 
perl). The PDF of the flux is also insensitive to the geometry 
of reionization. In addition, the PDF is heavily contaminated 
by noise and distorted by resolution. 

Dark gaps are defined to be continuous regions with 
optical depth above a threshold in spectra. Intuitively, the 
statistics of dark gap should contain the same information 
as leaks. However, the size of dark gaps are sensitive to 
the instrumental resolution, because higher resolution data 
contain more small leaks (e.g., Paschos & Norman 2005; 
Paperl), and they also are contaminated by observational 
noise. Moreover, dark gaps are from saturated absorptions, 
they are featureless and contains generally less information 
of non-saturated absorption. 

In this paper we made a statistical approach to 
Lya leaks. The purpose is to show that the statistical fea- 
tures of LyQ leaks would be effective tool to reveal the ori- 
gin of Lya leaks, and to probe the evolution of reionization. 
Similar to Lya absorption lines, Lya leaks have a rich set 
of statistical properties, such as the width of leak profile. 
Unlike dark gaps, the properties of Lya leaks can be defined 
through integrated quantities, which are less contaminated 
by resolution and noise. Moreover, Lya leaks are from re- 
gions of non-saturated absorptions and encode more infor- 
mation of reionization; therefore, the Lya statistics would 
provide more underlying physics of reionization than all the 
above-mentioned statistics. 

The paper is organized as follows. §2 describes the 
method to produce Lya absorption samples. §3 presents the 
statistical properties of Lya leaks with a uniform ionizing 
background. §4 analyzes the effect of inhomogeneous ioniz- 
ing background. Conclusion and discussion are given in §5. 



2 SIMULATION SAMPLES OF HIGH 

REDSHIFT LYa ABSORPTION SPECTRUM 

2.1 Method 

We simulate Lyman scries absorption spectra of QSOs be- 
tween z = 3.5 and 6.5 using the same lognormal method as 
those for low redshifts a ~ 2 — 3 (e.g., Bi et al. 1995; Bi 
& Davidsen 1997). In this model, the density field p(x) of 
the IGM is given by an exponential mapping of the linear 
density field 6o (x) as 

p(x) = po exp[5o(x) - o-o/2], (1) 

where ctq = (Sq) is the variance of the linear density field on 
scale of the Jeans length. Obviously, the 1-point PDF of p(x) 
is lognormal. In this model, the velocity field of baryon gas 
is produced by considering the statistical relation between 
density and velocity field (Bi & Davidsen, 1997; Choudhury 
et al. 2001; Veil et al. 2002). 

The dynamical bases of the lognormal model have grad- 
ually been settled in recent years. First, although the evolu- 
tion of cosmic baryon fiuid is governed by the Naiver-Stokes 
equation, the dynamics of growth modes of the fiuid can 
be sketched by a stochastic force driven Burgers' equation 
(Berera & Fang 1994). On the other hand, the lognormal 
field is found to be a good approximation of the solution of 
the Burgers' equation (Jones 1999). The one-point distribu- 
tion of the cosmic density and velocity fields on nonlinear 
regime are consistent with lognormal distribution (e.g. Yang 
et al. 2001, Pando et al. 2002). Especially, it has been shown 
recently that the velocity and density fields of the baryon 
matter of the standard ACDM model is well described by 
the so-called She-Leveque's universal scaling formula, which 
is given by a hierarchical process with log-Poisson probabil- 
ity distribution (He et al. 2006, Liu & Fang 2007). 

The simulation is performed in the concordance ACDM 
cosmological model with parameters flm = 0.27, h = 0.71, 
(78 = 0.84, and fib = 0.044. The thermodynamic evolution 
in the IGM is actually a rather complex process, because 
the nonlinear clustering leads to a multi-phased IGM. As 
shown in cosmological hydrodynamic simulations (e.g.. He 
et al. 2004), for a given mass density, the temperature of the 
IGM could have large scatters with differences up to two or- 
ders. Nevertheless, the equation of state in Lya clouds can 
be well approximated by a polytropic relation with 7 = 4/3 
( Hui & Gnedin; He et al. 2004). The neutral fraction /hi is 
obtained by solving the photoionization equilibrium equa- 
tion. The photoionization rate Fhi will be given in §2.2. We 
then construct synthetic absorption spectra by convoluting 
the neutral hydrogen density field with Voigt profiles. For 
each given z, the size of the simulation samples is Az = 0.3, 
and there are 2^* pixels in each simulation box. 



2.2 Redshift-Dependence of Photoionization Rate 

If the distribution of the IGM is uniform and the UV ionizing 
background is independent of redshift, the mean GP optical 
depth of Lya absorption should approximately increase with 
redshift as (l + z)''"''. The observations of dark gaps directly 
show that the GP optical depth undergoes a stronger evolu- 
tion at 2 ~ 6, and consequently, the UV ionizing background 
would decrease rapidly with redshift at 2 ~ 6. The strong 
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Figure 1. Redshift evolution of effective optical depth and neutral hydrogen fraction /hi with the photoionization rate given by eq.(2). 
The amplitude is taJien to be Fq = 5, 6 and 7. The data points are taJien from the Lya observation of Fan et al. (2006). 



evolution scenario of the UV ionizing background is sup- 
ported by a number of simulations or semi-analytical models 
of reionization (e.g., Raaoumov et al. 2002; Gnedin 2004; Oh 
& Purlanetto 2005; Pascho & Norman 2005; Wyithe & Loeb 
2005; Kohler et al. 2007; Gallerani et al. 2006). It has been 
found that an evolution of photoionization rate as follows 
can fit the strong redshift evolution of the GP optical depth 
(Paper I) : 

Vuiiz) = Fo exp{-[(l + z)/{l + 3.2)f -^}, (2) 

which is in units of 10~^^ s~^. Note that the power index 
2.4 in equation (2) is little different from the one used in 
Paper I because we use a different To in this paper (also see 
below) . 

With cq.(2), wc calculate the redshift dependencies of 
neutral hydrogen fraction /hi and the effective optical depth, 
Teff = — ln(_F), where F is the mean of transmitted flux. 
The results are plotted in Figure 1. The data points for Teff 
are taken from Lya observations of Fan et al. (2006). For 
best fitting, the amplitude Fo is in the range 5-7, which can 
be considered as the allowed range of Tm{z). In this paper, 
we will use Fo = 6 as the fiducial photoionization rate. It is 
interesting to note that fm approaches to ~ 0.1 at redshift 
z ~ 10, which is consistent with the electron scattering opti- 
cal depth given by the data of CMB polarization of WMAP 
III (Page et al. 2007). 

It should be pointed out that the assumption of To = 2 x 
lO'' K (§2.1) is well reasonable at z < 5 (e.g. Hui & Gnedin 
1997; He et al 2004) and may still be applicable at 2: ~ 6 if 
the mass averaged neutral fraction of hydrogen is not larger 
than 10~^, and the photon heating rate is small. However, 
at higher redshift, say z > 6, the temperature To might be 
redshift-dopcndcnt. Yet, no proper information on To{z) is 
available at high redshift, and this leads to uncertainty of 
the model. Fortunately, in photoionization equilibrium, the 
neutral fraction /hi depends mainly on a degenerate factor 
rHi(2)To '^^(z). Thus, the problem with the uncertainty of 
To{z) can be overcame if wo use the combined parameter of 
rHi(2)[To(«)/2 X 10"]°-^^ to fit the data. In the range z < 6, 
this parameter actually is Fhi(z); in the range of 2: > 6, it 
is different fi-om Fhi(^;) by a factor of [To(^;)/2 x 10''1°-''^ 



Thus, the rodshift-ovolution of Fhi(z) would bo slower than 
eq.(2) if To{z) is less than 2 x 10** K at higher redshifts. 



2.3 An Example of Lya Absorption Spectrum at 

z = 6 

As an example of Lya leaks, we plot a simulated sample of 
Lya absorption spectrum at z = 6 with a uniform ionizing 
background in Figure 2, which shows the transmitted flux 
F, the density p of baryon gas, the bulk velocity v, and 
the neutral hydrogen fraction /hi. As expected, the mean of 
transmitted flux is very small, about 0.004, and corresponds 
to an effective optical depth 5.5. Nevertheless, we see spiky 
features with the transmitted flux F as large as 0.15. They 
are leaks. 

At low redshifts z < 5, the Lya forests in QSO's spectra 
have a spectral filling factor significantly less than one and 
can be decomposed into individual Lya absorption lines. At 
redshifts 2 > 5, it is meaningless to decompose the spectra 
into individual lines since almost the whole spectra are ab- 
sorbed completely. We note, however, the transmitted leaks 
look like emission features upon the dark background, and 
the absorption spectra can be decomposed into individual 
"emission lines", i.e., Lya leaks. 

Comparing the top, the second and bottom panels, we 
see that all the leaks comes from the regions with mass 
density less than 0.3 of the mean mass density of baryon 
gas. The neutral fraction for leaks is /hi ~ 2 x 10~^, which 
yields a GP optical depth ~ 2.5 for overdensity 0.25 and a 
F ~ O.l, while the mean neutral fraction of the entire exam- 
ple is about 7 x 10~^, which is high enough to produce dark 
gaps. The column density of neutral hydrogen of the leaking 
features is mainly in the range of 10^^-10" cm , which are 
the non-saturated absorption regions, and therefore, leaks 
can come from regions where no enough neutral hydrogen 
to produce complete absorptions. 

Similar to very high density clouds, p/p ^ 1, the re- 
gions with very low density p/p <C 1 arc rare events in the 
cosmic clustering. Therefore, leaks may provide valuable test 
on models of clustering. For instance, the lognormal distri- 
butions are long tailed in both high and low density sides, 
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Figure 2. An example of simulated spectrum at z = 6. It shows the transmitted flux F, the density p of baryon gas, the velocity v, and 
the fraction of neutral hydrogen /hi from top to bottom, (p) is the mean of baryon matter. 



and it contains more high density events as well as low den- 
sity events than Gaussian model. It should be emphasized 
that, once the photoionization rate Fhi is determined from 
the GP optical depth, the statistical property of the sample 
shown in Figure 2 doesn't contain free parameters. These 
samples have been successfully used to explain the follow- 
ing observations: 1.) the large dispersion of the GP optical 
depth; 2.)the PDF of the flux, and 3.) the evolution of the 
size of dark gaps (Paper I). Now we use them to study the 
statistical properties of Lya leaks. 



3 STATISTICAL PROPERTIES OF LYa LEAKS 

The Lya absorption lines at low redshifts are described by 
well defined quantities, such as the equivalent width, FWHM 
(full width half maximum) , and the Voigt profile, all of which 
are easily related to physical interpretations. Since the trans- 
mitted leaks look like emission features, one can decompose 
the spectra into individual "emission lines" and describes 
it by quantities similar to the Lya absorption lines. In this 
section we study the statistical distribution of Lya leaks in 
the model with a uniform ionizing background. The effects 
of inhomogeneous ionizing background will be discussed in 
next section. 



0.100 
W(A) 



Figure 3. Relation between W and Wh of Lya leaJcs at z = 6. 



3.1 Profile of Leaks 

The center of a leaJj is identified as the maximum of trans- 
mitted fiux, and the boundaries are two nearby positions 
around the center where the flux falls down to zero or to 
noise level (we take it to be F = 0.001 in this paper). The 
properties of leaJis can be measured by two quantities: an 
equivalent width W, which is the total area under the pro- 
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file of a leak, and a half width Wh, which is defined as the 
width around center within which the area under the profile 
of the leak is equal to haJf of the total area of the leak. 

The equivalent width W and the half width Wh are 
two independent measurements of Lya leaks. This point can 
be seen in Figure 3, which gives the value of W and Wh 
for each leak at z = 6. It does not indicate a significant 
correlation between W and Wh, especially in the region of 
> 0.1 A and Wh > 0.1 A. Figure 3 also shows that the 
distribution of Wh has a lower limit 0.2 A, which is due to 
the Jeans scale used for smoothing the sample. On the other 
hand, the equivalent width W distributes in the range from 
0.001 A to 5 A. Obviously, one is unable to introduce two 
independent quantities for characterizing dark gaps. 

The profiles of 100 randomly sampled leaks with W = 
0.01, 0.1, and 1 A are displayed in Fig.4. The tails of the pro- 
files in the three panels of Fig.4 look like the Lorentz profile, 
and of course, they do not have the meaning of the natural 
width of an absorption line. The profiles of leaks for a given 
equivalent width W have very large dispersions; for exam- 
ple, for W=0.lA, the fiux covers a range from F = 0.03 to 
0.15, and Wh can be 0.4 to 4 A. As the leaks are formed out 
of the two neighboring complete absorption troughs, which 
depend on inhomogeneities of density, velocity, and temper- 
ature fields, the large scatter of the profiles is expected. 

Figure 4 also shows that some leaks may have multi- 
ple local maximums above the noise level. For clarity and 
easy-operating, we treat them as one leak. The current ob- 
servational resolution is of the order of lOkm/s, which corre- 
sponds to AA ~ 0.28A in observer's frame, or about 13 pixels 
of simulation. Therefore, the observed Lya leaks would have 
a resolvable width. 

3.2 Number Density Distributions of Lya Leaks 

Similar to Lya forests, we can define the cumulative num- 
ber densities n(> W,z) and n(> Wh,z) of leaks as the 
number of leaks with widths larger than a given W and 
Wh at z per unit z, respectively. The differential number 
densities are n{W,z) = dn{> W,z)/dW and n{WH,z) = 
dn{> WH,z)/dWH- It should be pointed out that statistics 
of W and Wh are not the same as the largest peak width 



statistics proposed by Gallerani et al. (2007), which consid- 
ered only the largest peak width. A peak may contains more 
than one leaJis, i.e. leaJi statistics describe the details of the 
leaking area. The mean transmitted flux at z within dz is 
F — J^n{W,z)WdW. We calculate the number densities 
of leaks in redshift range z = 5 — 6.5. In each redshift re- 
gion we produce 100 light-of-sight samples to calculate the 
density functions. The results are shown in Figure 5. The 
errors are estimated by Jackknife method, i.e., the variance 
over 5 subsamples, each of which contains 20 light-of-sight 
samples. 

The number density n{W, z) shown in the up-left panel 

of Figure 5 are similar to a Schechter function: they follow 
a power law at small W and have a cut-off at large W. The 
distribution of n{WH, z) at Wh < 0.5A declines with de- 
creasing Wh, this is because of the Jeans length smoothing. 

The slope of n{W,z) and n{WH, z) are smaller for small 
redshift. It means the lack of low density voids with small 
size. That is, the increase of voids of small size is less than 
voids of large size. This trend can also be seen from the flat- 
tening of the cumulative density distributions n(> Wh,z) 
and n(> W, z). 

The redshift-evolution of the rmmber densities n,(> 
W, z) and n(> Wh, z) of leaks for W, Wh =0.5, 1.0, and 3.0 
A are shown in Figure 6. As has been seen in Figure 5, the 
number densities dramatically decrease at higher redshifts. 
The evolution is more rapidly for large leaks: the number 
density n(> W = 0.5A, z) drops by a factor of ~ 5 when 
redshift increasing from 5 to 5.8, while n(> W = 3A, z) 
drops by a factor of ~ 60. The evolution trend is the same 
for the number density n{WH, z). From the error bars of Fig- 
ures 5 and 6 we see that the predicted features of leaks would 
be able used to compare with observed data set containing 
20 or more light-of-sight samples with the similar quality as 
simulation. 

3.3 Effects of Resolution and Noise 

In this section, we study the observational and instrumental 
effects on the statistics of leaks. Since both W and Wh are 
defined through the area under the profile of leaks, the ef- 
fects of resolution and noise would be small. To simulate the 
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Figure 6. Rcdshift evolution of number densities n(> W,z), n(> W,z), n(Wi{ , z), and n(> Wu^z) of leaks with W, Wu =0.5, 1.0, 
and 3.0 A. The error bars are the l-u range given by Jackknife estimation, in which each subsample contains 20 lines of the absorption 
spectrum. 
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Figure 7. Effects of resolution and noise on number densities n(> 
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Figure 8. Number densities n(> W, z), and n(> Wh, z) for sample 
are the same as Figure 5. 

observational effects, we bin the original data to a coarse grid 
corresponding to a resolution of 20000, and we add Gaussian 
noises with signal-to- noise ratio S/N=3 on binned pixels. 
The number densities of leaks for the noisy binned samples 
are shown in Figures 7. 

The effect of binning and noise is very small for W: 
the original plot of W actually is the same as the plot of 
binned W, and the distribution of noisy W is affected only 
when W < O.OlA. The effect of binning for Wn is also small 
on scales larger than the binning length. As expected, the 
noise effects for Wh arc significant for Wn < 0.4 A. The 
noise effects for Wh arc even smaller if we smooth the noisy 
sample. This is very different from the PDF of the flux and 
the size of dark gaps, both of which are heavily contaminated 
by instrumental resolution and observational noises. 

One can compare the uncertainty of Fq with the effect 
of noises. Figures 8 shows the number densities for the UV 
background amplitude Fo = 5, 6, and 7 (cq.2), which rep- 
resent the allowed range of To. Different from noises, the 
difference of Fo will cause uncertainty in the whole ranges 
of W and Wh- The uncertainties of number densities are 
within a factor of 2 when the amplitude Fo changes from 




Wh(A) 

\z) (left panel) and n{> Wh,z) (right panel) for leaks at 2 = 6. 




at 2 = 6 with Fo = 5, 6, 7 xlO s ^, respectively. The error bars 



7 to 5. These uncertainties essentially are from the mass 
density perturbations with long wavelengths (Paper I). The 
error bars from the scattering of light-of-sight samples is 
also shown in Figure 8. Therefore, the scattering of > 20 
light-of-sight samples actually is less than the uncertainty 
of Fo. 



4 Lya LEAKS OF INHOMOGENEOUS 
IONIZING BACKGROUND 

In the early stage of reionization, ionizing photons are 
mainly in ionized patches around UV ionizing sources, and 
therefore, the spatial distribution of ionizing background 
is highly inhomogeneous and has patchy structures. When 
the ionized patches spread over the whole space, the ion- 
izing background become uniform or quasi-uniform, and so 
the ionizing background underwent an inhomogcncous-to- 
uniform transition during the reionization. In this section 
we study the effects of inhomogeneous ionizing background 
on Lya leak statistics. 
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4.1 Models of Inhomogeneous Ionizing 
Background 

The first question is how to model the inhomogeneous spa- 
tial distribution of the ionizing photon field. To our prob- 
lem, the most important property is the relation between 
the fields of mass density and ionizing photon background 
or photoionization rate. If the spatial fluctuation of pho- 
toionization rate rHi(x,2;) is statistically uncorrclated with 
the density field p(x, «), the reionization of IGM will sta- 
tistically be the same as a uniform ionizing background, re- 
gardless of the details of the fluctuation of rHi(x, z). In this 
case, the only effect of the fluctuations would be to yield a 
larger variance in relevant statistics. However, as shown in 
last sections, the uncertainty of the leak statistics is already 
large even when the ionizing background is uniform, and 
therefore, one would not be able to distinguish the fluctu- 
ating photon field from inhomogeneous density field if they 
are uncorrelated. 

Although many simulations on the UV photon field at 
the epoch of reionization have been done, there still lack the 
results of the correlation between photon and density fields. 
In this context, we will consider the following four models 
on the statistical relation between the inhomogeneous fields 
of photon and density, which arc mainly based on physical 
consideration of different stage of reionization. 

Model 0. The photoionization rate is spatially uniform. 
It corresponds to the post-overlapping stage of reionization. 
This model has been used in last two sections. 

Model I. The photoionization rate at a give redshift is 
assumed to be proportional to the density field of IGM, 
Fhi = Tip, Fi being a constant. This model is motivated by 
the so-called inside-out scenario: high density regions around 
UV sources are ionized first, and is most probable at the 
early stage of reionization. 

Model II. Just opposite to model I, the photoionization 
rate at a give redshift is assumed to be inversely propor- 
tional to the density field of IGM, Fhi = Fnp^"*^, Fn being 
a constant. This model comes from the so-called outside- 
in scenario: under-dense regions are ionized first (Miralda- 
Escude et al. 2000), which is applicable at the late stage of 
reionization (Purlanetto & Oh 2005). 

In order to fit the observed effective optical depth r^ff 
at redshift 2: = 6, we take the following parameters: Fi = 
Fo X 3.53, Fii — Fo/2.77 where Fo is the photoionization 
rate at z = 6 for uniform ionizing background [eq.(2)]. 

Model III or patch model. This model corresponds to 
the stage before the overlapping of ionized patches. The 
IGM is almost fully neutral except for some isolated ionizing 
patches. We model the ionized patch as a Stromgren sphere: 
the neutral hydrogen fraction is equal to within the patch 
and equal to 1 outside the patch. The mean radius of the 
Stromgren sphere is assumed to be Rs — 1.8 comoving Mpc, 
which corresponding to an UV photon source with luminos- 
ity L = 5 x 10*^ erg s~^ with active time 10^ year and a 
spectrum. To fit the observed optical depth, we found there 
should be 2 patches every simulated box {Az — 0.3). This 
number is actually consistent with the following estimation 

dN _ ■kR14>{L)c 

dz H{z) ' ^ ' 

where (jyi^) the comoving luminosity function, i.e. the 3-D 



number density of sources with luminosity L. Using 0(5 x 
10*^) ~ 1.5 X 10"^ Mpc"'' (Bouwens et al. 2006), we have 
dNIdz ~ 7 and {dN/dz)Az ~ 2. 

4.2 Lya Leaks of Inhomogeneous Ionizing 
Background 

For the four models described in the last section, we recalcu- 
late the transmitted flux with the same underlying density 
and velocity fields of fig2. The results are shown in Figure 
9. As mentioned above, the mean of transmitted flux for all 
models takes approximately the same value of 0.004, which 
corresponds to an effective optical depth ~ 5.5. However, 
Figure 9 shows clearly different behaviors of Lya leaks in 
various models. 

For models I and II, the leaks appear exactly at the 
same positions as model except some small leaks, which 
are more prominent for model I. In other words, all the leaks 
of models 0, I and II are from lowest density voids. Therefore, 
the three models have the same distribution of dark gaps, 
and it is impossible to discriminate among these models with 
the dark gap statistics. 

However, the profiles of the leaks of model 0, 1 and II are 
statistically different from each other. For model I, the leaks 
generally have larger width and lower height than model 0, 
while for model II, the width of leaks generally is narrower 
than model 0, and the height of leaks is larger than model 
0. The reason is straightforward. Comparing with model 0, 
model I gives a higher /hi at low density and lower /hi 
at high density. This leads to lower amplitude and broader 
width. For model II, the effect is just opposite to model I 
and yields higher amplitude and narrower width. 

In the patch model, all the leaks come from ionized 
patches, within which ionizing sources arc enclosed, and so 
the internal information of ionized patches can be inferred 
from the leak statistics. The size of dark gaps is actually 
given by the distance between ionized patches. Generally, 
Lya leaks in the patch model have a maximum flux ~ 0.6, 
which is higher than the maximum flux of other models, 0.2. 
This is because the neutral fraction within the Stromgren 
sphere is much less that than other models. This behavior 
is similar to the so-called proximity effect of QSOs at low 
redshift (e.g., Ranch 1998). Thus, the statistics of the max- 
imum flux of Lya leaks may be used to reveal the patchy 
origin of Lya leaks. 

As we sample the size of patches along the sight of light 
according to the impact probability, the size of patches is 
smaller than Rs. If the size is too small, the ionized patch 
will be opaque to Lya photons due to the damping wing 
of the surrounding neutral hydrogen absorption (Miralda- 
Escude 1998). Obviously, it explains why only one leak is 
apparent in the patch model as displayed in Figure 9. 

4.3 Number Densities of Lya Leaks of 
Inhomogeneous Ionizing Background 

We now calculate the number densities n{W,z), n{> W,z), 
n{Wu,z) and n(> Wh,z) for three inhomogeneous ionizing 
background models at a = 6, and the results are shown in 
Figures 10. We can see from Figure 10 that the effects of 
different models on the number densities of W and Wh are 
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Figure 9. Transmitted flux F given by uniform ionizing background (top panel), the Model I (second panel), the Model II (third panel) 
and the Model III (fourth panel). The density field p of baryon gas (bottom panel) is the same as Figure 2. 



different. That is, although all the models of 0, I, II and 
III give the same mean effective optical depth, their leak 
distributions are different. For clarifying, in Figures 10 we 
show only error bars for the curves of model 0. It would be 
enough to show that these results can be tested with data 
set containing about 20 light-of-sight data. 

First, for model I, both of the number densities n{W, z) 
and n(> W, z) have only small deviations from model 
0. However, n{WH,z) are significantly different from their 
counterparts of model 0. The number density of leaks with 
Wh > lA is much more than model 0. The basic feature of 
model I is to increase the number of leaks with large Wu as 
shown in Figure 9. 

Second, for model II, we see once again that the number 
densities n{W,z) and n(> W,z) have only small deviations 
from model 0. Yet, the number densities n{Wu,z) and n{> 
Wii,z) of the model II are systematically lower than model 
0. Therefore, the basic feature of model II is to keep the 
total area under the profile of leaJis almost unchanged, but 
the widths of leaks are significantly narrowed. 

Finally, the behavior of patch model is very different 
from models 0, I and II. As expected, the patch model 
yields more leaks with large W, and n{W, z) shows a bump 
around 14-' = 3 — 5 A, which characterize the area of the 
ionized patches. On the other hand, the number densities of 
n(WH,«) and n(> Wyi,z) are lower than model because 
the characteristic size of ionized patches is less than that 



of low density voids. It is interesting to note that for the 
patch model, the tails of n{W,z), n(> W,z) and n{Wn,z), 
n{> Wn,z) are quite different from each other. Generally, 
the tails of n{W, z) and n{> W, z) can extend to as large as 
W ~ SA, while for n{Wii,z) and n(> Wh,z), there are no 
tails higher than 3 A. It results partially from the damping 
wing effect (Miralda-Escude 1998). 

In summary, the statistical properties of leaks with re- 
spective to W and Wh are sensitive to the details of ionizing 
photon field. Combining the distribution of W and Wh, the 
Lya leaks would be able to probe the origin of themselves, 
and thus reveal the ionization state of IGM, the iuhomo- 
geneity of ionizing background, and the evolution stage of 
reionization. 

It should be pointed out that we considered only the 
patches of the HI regions around high redshift galax ies. The 
HI regions around quasars or the proximity effect would also 
be the patches of leaking. The mean luminosity of quasars 
probably is higher than galaxies, and therefore, the above- 
mentioned feature of n{WH,z) and n(> Wh,z) would be 
more prominent for the patches of quasars. The damped 
Lya absorption system is important for modeling low red- 
shift Lya absorption. Since these systems have high column 
density of neutral hydrogen, it will not contribute to leaking 
either in density void models or patch model. 
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Figure 10. Number densities n{W,z),n{> W,z), n{WH,z), and n(> Wh,z) of leaks for tliree fluctuating ionizing background models 
at z = 6. The error bars of the uniform background model are taken from Figure 5. 



5 DISCUSSION AND CONCLUSION 

We show that the Lya absorption spectra of UV photon 
emitters at redshifts around a ~ 6 can be decomposed into 
Lya leaks, which come from cither lowest density voids or 
ionized patches containing ionizing sources. The Lya leaks 
are well defined and described by the equivalent width W 
and the width of half area Wh- Since both W and Wh are 
defined through integrated quantities, the distributions of 
Lya leaks in terms of W and Wh are stable with respect to 
observational noises and instrumental resolution. Although 
the number densities n{W, z), n(> W, z), n{WH, z), and n(> 
Wh, z) evolve very rapidly at redshift z ~ 6, these statistics 
are measurable up to ^ = 6.5. 

If the Lya leaks come from lowest density voids, the 
IGM should be still highly ionized and the ionizing back- 
ground is almost uniform; in contrast, if the leaks come from 
isolated ionized patches, the ionizing background should bo 
inhomogeneous, and the reionization is still in the overlap- 
ping stage. Therefore, the origin of Lya leaks is crucial to 
understand the history of reionization. 

We show that the statistics with W and Wh are sensi- 
tive to the origin of Lya leaks, because the Lya leaks are sen- 
sitive to the correlation between photon and density fields. 
Based on physical consideration of reionization, we stud- 
ied four phenomenological models of the photon field: the 



uniform ionizing background (model 0) , the photoionization 
rate Fhi is proportional to the density p (model I), Fhi and p 
are anti-corrclated (mode II), and ionized regions only given 
by Stromgren sphere around ionizing photo sources (patch 
model). We found that, although all the four models can fit 
the observed mean and variance of optical depth at ^; ~ 6, 
the width distribution of Lya leaks show very different be- 
haviors. 

For model 0, 1, and II, most of Lya leaks are from lowest 
density voids, and the distribution of dark gaps are similar. 
However, the properties of individual Lya leaks are differ- 
ent. Model I gives broader width Wh than model 0; model 
II gives narrower width Wh than model 0. For patch model, 
the Lya leaks are mainly from ionized sphere, and they gen- 
erally have higher maximum of transmitted flux than other 
models. There is a bump in the distribution of the equivalent 
width W , which characterize the intensity of UV photons 
within ionized patches. Lya leaks from ionized patches will 
provide the information of ionized patches, such as their in- 
tensity and size, and constrain the properties of UV sources 
that contribute most of ionizing photons of the reionization. 

Finally, we point out that similar analysis is also appli- 
cable to the reionization of Hell. The optical depth of Hell 
Lya evolves rapidly at « ~ 2 — 3. It reaches ~ 5 at « ~ 3. 
That is, the evolution of Hell reionization at « ~ 3 would 
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be similar to that of H at a ~ 6. Observed sample of HoII 
Lya absorption indeed show structures similar to the leaks 
of hydrogen absorption spectrum (Zheng et al. 2004). One 
can also draw the information of background photons at the 
energy band of Hell ionization. 
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